Lyapunov Generation of Entanglement and the Correspondence Principle 
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We show how a classically vanishing interaction generates entanglement between two initially 
nonentangled particles, without affecting their classical dynamics. For chaotic dynamics, the rate 
of entanglement is shown to saturate at the Lyapunov exponent of the classical dynamics as the 
interaction strength increases. In the saturation regime, the one-particle Wigner function follows 
classical dynamics better and better as one goes deeper and deeper in the semiclassical limit. This 
demonstrates that quantum-classical correspondence at the microscopic level requires neither high 
temperatures, nor coupling to a large number of external degrees of freedom. 
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In the decades since its inception, no observed phe- 
nomenon, nor experimental result ever contradicted 
quantum theory. Yet, the world surrounding us, though 
being made out of quantum mechanical building blocks, 
behaves classically most of the time. This suggests that, 
one way or another, classical physics emerges out of quan- 
tum mechanics. Today's common understanding of this 
quantum-classical correspondence is based on the real- 
ization that no finite-sized system is ever fully isolated. 
It is then hoped that a large regime of parameters exists 
where the coupling of the system to external degrees of 
freedom (to be called the environment from now on) de- 
stroys quantum interferences without modifying the sys- 
tem's classical dynamics. Indeed, such a coupling usually 
induces loss of coherence on a time scale much shorter 
than it relaxes the system P, Q . 

The standard approach to decoherence starts from a 
master equation valid in the regime of weak system- 
environment coupling 0, 0] . The master equation deter- 
mines the time-evolution of the system's Wigner function 
M^(p,q) = {2nh)-'^J dxexp[ipx]p(q + x/2, q - x/2) {p 
is the system's density matrix) as 

dtW = (h,w] + Y^ ''^^ a2"+V92"+iT/F 

* I ' / 22"(2n-f 1)! « P 

n>l ^ ' 

+2-idp{pW) + DdlW. (1) 

The first term on the right-hand side of Eq. ^ is the 
classical Poisson bracket. The second term, written here 
for the case of a momentum- independent potential V{c\), 
exists already in closed systems and generates quantum 
corrections to the dynamical evolution of W . This term 
starts to become comparable to the Poisson bracket at 
the Ehrenfest time te — A~"'^lnA^, where A is the Lya- 
punov exponent of the classical dynamics, and N the 
size of the system's Hilbert space. The two terms on 
the second line of Eq. are induced by the coupling to 
the environment. In the limit of weak coupling, 7-^0, 
but finite diffusion constant, D oc — Cst, the fric- 
tion term vanishes, leaving the classical dynamics unaf- 
fected. This requires to consider the high temperature 
limit. Simultaneously, for large enough D, the momen- 
tum diffusion term induces enough noise so as to kill the 



quantum corrections before they become important. The 
time-evolution of W is then solely governed by the clas- 
sical Poisson bracket, that is to say, classical dynamics 
emerges out of quantum mechanics. Refs. .Jj dj] provided 
for a numerical illustration of this scenario. 

Our motivations in this article are as follows. First, it 
is unclear how generic the above scenario is, since it is 
based on a master equation derived under restrictive as- 
sumptions, for instance on the environment, the dimen- 
sionality of the system or the strength of the coupling 
between system and environment Also, it formally 

requires to consider infinite temperatures. Moreover, and 
with the specific exception of the kicked harmonic oscil- 
lator investigated in Ref. 3, there is not much analyti- 
cal understanding of the decoherence process in generic 
dynamical systems, i.e. except for the regular case, mas- 
ter equations are usually integrated numerically. Second, 
claims have been made of an entropy production due to 
the coupling to environmental degrees of freedom gov- 
erned by the system's exponent A 0,13, without clear 
analytical derivation, nor strong numerical evidence |^. 
A Lyapunov decay of the fidelity has recently been ana- 
lytically predicted ff| and numerically observed , how- 
ever, decoherence and fidelity are two different things, 
especially in the generic situation where the system and 
environment Hamiltonians do not commute with the cou- 
pling Hamiltonian 9, 10]. 

We revisit these issues and consider two interacting 
quantized dynamical systems. Entanglement generation 
between two particles has already been considered in 
Refs. US m III El El- All results to date are con- 
sistent with the scenario proposed in Ref. Ell , according 
to which bipartite entanglement results from two contri- 
butions: (i) a quantum-mechanical one, which depends 
on the coupling strength between the two systems, and 
(ii) a dynamical one, which, in chaotic systems, is deter- 
mined by the total system's spectrum of Lyapunov expo- 
nents. The entanglement generation rate is given by the 
weakest of the coupling strength and the Lyapunov expo- 
nent. It has to be pointed out that this picture holds in 
the regime of classically weak but quantum-mechanically 
strong coupling (this will be made quantitative below). 
For regular systems, entanglement generation is slower 
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than for chaotic ones, typically power-law in time |lCtll4| . 

The purpose of this letter is threefold. First, we ad- 
dress the problem of decoherence and bipartite entan- 
glement from a microscopic point of view, i.e. without 
relying on effective differential equations. This allows 
for a clear identification of the regime of validity of our 
theory. Second, we give strong numerical evidences for 
the existence of a Lyapunov regime of entanglement in 
bipartite systems (the numerical evidences presented in 
Ref. ^3 were challenged in Ref. Third, we dis- 

cuss our results from the point of view of the quantum- 
classical correspondence, and present numerical phase- 
space dynamics results showing that this correspondence 
is fully achieved in the regime of Lyapunov entanglement. 
This is, we believe, the first clear microscopic illustra- 
tion of the quantum-classical correspondence in a generic 
chaotic system. 

As our starting point, we consider the Hamiltonian 

n = Hi®l2^h®H2 + hU. (2) 

We take chaotic one-particle Hamiltonians Hi_2- We 
specify that the interaction potential hi is smooth, vary- 
ing over a distance much larger than the particles' de 
Broglie wavelength cr, and that it depends only on the 
distance between the particles. Planck's constant in front 
of U in Eq. Q makes it explicit that we consider a 
semiclassically vanishing two-particle interaction, i.e. the 
classical Hamiltonian corresponding to Eq. Q does not 
couple the two particles. Our goal is to calculate the 
purity V{t) = Tri[p^(i)] of the reduced density matrix 
pi(t) = Tr2p(t). In the situation we consider of a unitary 
two-particle dynamics acting on an initially pure two- 
particle state, V^t) is a good measure of entanglement, 
which varies between 1 for product states and for max- 
imally entangled states. The calculation proceeds along 
the lines of Ref. and here we only sketch it. A similar 
semiclassical approach has been applied to a stochastic 
Schrodinger equation in Ref. ■ 

In the initial two-particle product state we take, 
each particle is in a Gaussian wavepacket tpi 2 (y) = 

(w2)-'i/4exp[zpi,2 • (y - v,a)/n - |y - ri,2|V2^']- 
The two-particle density matrix evolves according to 
pit) — exp[— iTit/fi,]!?/;!; ?/'2}('0i; "021 exp[i7it/fi] This time- 
evolution is evaluated semiclassically by means of the 
semiclassical two-particle propagator 

(xi,X2|e-«*/''|yi,y2) ^ C,\/?e^[(^=+^=')/^+5„.,] 

s,s' 

(3) 

which is expressed as a sum over pairs of classical trajec- 
tories, labelled s and s' , respectively connecting yi to xi 
and y2 to X2 in the time t. Because of our assumption of 
a semiclassically vanishing coupling, these classical tra- 
jectories are determined by the one-particle Hamiltoni- 
ans. Each pair of such trajectories gives a contribution 
weighted by Cs,s', the inverse of the determinant of the 
stability matrix on s and s', and oscillating with one- 
particle (denoted by and Ss') and two-particle (de- 



noted by Ss,s' = /odi'iY[qi^^(t'),qi'^ (*')]) action inte- 
grals accumulated by the first and second particles along 
s and s' respectively. In the regime of semiclassically 
vanishing coupling we consider, the onc-particle actions 
generate much faster oscillations than their two-particle 
counterpart. Accordingly, our approach relies on station- 
ary phase conditions imposed on the one-particle actions. 
In Eq. (PJ, Maslov indices have been omitted since they 
drop out of the calculation. 

To leading order in Hcs — 2tt/Ni,2 {Ni is the size of the 
system's Hilbert space), our semiclassical calculation 
gives the time-evolution of the purity as 

P{t) ~ ^ a,e(t > r,)exp[-Aji] +exp[-2rt] 

4=1,2 

+e{t>Ti'^)N^' + eit>4'^)N^\ (4) 

The first, classical term decays with the Lya- 
punov exponents Ai^2 ^l^l- It does not exist at 
short times, t < Ti — X^^ \n[\i/a^Gi], {Gi — 

/dt'(a«W[qi^)(0),qi?)(0)] a«W[qi^)(f),qi?)(t')])), and 
has prefactors — 0{1). The second term is the stan- 
dard, interaction-dependent quantum term with F = 
j;;dt'{U[qi'\0),q^^\0)] U[q^^\t'),q^^\t')]), assuming a 
fast decay of correlations. Being given by the correla- 
tor of a classical potential evaluated along classical tra- 
jectories, F does not depend on h. Finally, the third 
and fourth, saturation terms in Eq. Q set in after 
r«=ArilniV,. 

The validity of our approach is given hy 62 < F < i?2 , 
where B2 and 62 = B2 / {N1N2) are the two-particle band- 
width and level spacing respectively 01 . In this range, U 
is quantum-mechanically strong as individual levels are 
broadened beyond their average spacing, but classically 
weak, as B2 is unaffected by W [3, 0- We note that 
our semiclassical approach preserves the properties of the 
density matrix Tri[pi(t)] = 1, pi = p|, as well as the 
symmetry TrM{t)] = Tr2[pl{t)]. 

Eq. expresses the decay of V{t) as a sum over dy- 
namical, purely classical contributions, and quantal ones, 
depending on the interaction strength. Because the de- 
caying terms are exponential and have prefactors of order 
unity, one has for i > , ri.2 

7'(t) 2±exp[-min(Ai,A2,2F)i]+Arfi+7V2-\ (5) 

Eq. lO reconciles the results of Refs. |0| and Its 
regime of validity, 82 — B2/ {N1N2) < F < i32, is para- 
metrically large in the semiclassical limit iVi.2 — > 00. The 
same approach also applies to regular systems, in which 
case the exponentially decaying Lyapunov terms are re- 
placed by power-law decaying terms [lol IT^ . 

We now discuss the connection of our main result, 
Eq. ijSJl, to Eq. The purity measures the weight of off- 
diagonal elements of pi(t), and hence of the importance 
of coherent effects. In the regime 2F Ai = A2, P(t) 
reaches its minimal value at the Ehrenfest time. Thus, 
quantum effects [the second term on the right-hand side 
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FIG. 1: Main plot: Purity of the reduced density matrix 
for A^i = iV2 = 512, Ki = K2 £ [4, 12], and e = 4 giving 
2V — 13.6 S> Ai — X2- Data have been calculated from 20 
different initial states. The time axis has been shifted by the 
onset time ri (see text) and rescaled with Ai G [0.5,1.35]. 
The full line indicates oc exp[— Aif], and the dashed line gives 
the asymptotic saturation T'(oo) = 2A'^j~^. Inset: Purity 
for Ki = K2 = 5.09 for e = 0.2 (circles), 0.4 (squares), 0.8 
(diamonds), 1.6,2,3 and 4 (triangles). 

of Eq. 1^] are killed before they have a chance to ap- 
pear. In that regime, one therefore expects the quantum- 
classical correspondence to become complete in the semi- 
classical limit Ni^2 — > 00. We will show below numerical 
evidences supporting this reasoning. 

To numerically check our results, we consider Hamil- 
tonian Q for two coupled kicked rotators [3 

H, - p'^/2 + K,cosix,) J2S{t-nT), (6a) 

n 

U = e sin(xi - X2 - 0.33) ^ S{t - nT). (6b) 

n 

The interaction potential U is long-ranged with a 
strength e and acts at the same time as the kicks. Upon 
increasing Ki the classical dynamics of the i^^ particle 
varies from fully integrable {Ki — 0) to fully chaotic 
[Ki > 7, with Lyapunov exponent Ai ~ \n{Ki/2)]. For 
1 < Ki < 7 the dynamics is mixed. We will vary 
Ki^2 € [3,12] to get a maximal variation of Ai, while 
making sure that both tpi and ^'2 he in the chaotic sea. 
We follow the usual quantization procedure on the torus 
x,p € (— 7r,7r). The bandwidth and level spacing are 
given by B2 ^ 4:Tt, 62 ~ 47r/iV^, and we numerically 
extracted T ~ 0.43e^ from exact diagonalization cal- 
culations of the local spectral density of states. The 
time evolved density matrix is computed by means of 
fast fourier transforms |l7j . The algorithm requires only 
0{N In N) operations, which allowed us to reach system 
sizes up to iVi_2 = 2048, more than one order of magni- 
tude larger than any previously investigated case. 

The behavior of V{t) is shown in Fig. First, it is 
seen that as e increases, the rate of entanglement gen- 
eration also increases, up to some value ec, after which 
it saturates. We have found that (i) prior to saturation. 



'P{t) decays exponentially with a rate « 0.85e^, provided 
r = 0.43e2 > S2 = Att/N'^ is satisfied, and that (ii) Cc 
behaves consistently with Eq. (jSJ. Second, Fig. 1 shows 
how V{t) behaves for fixed e > Cc upon variation of the 
Lyapunov exponents Ai = A2. The rescaling of the time 
axis t — > Alt allows to bring together six curves with 
Ai G [0.5, 1.35], varying by almost a factor three. Third, 
Fig. 1 shows that in the chaotic regime considered here 
with Ni = N2, V{t -> cx)) = 2N^^. These numerical 
data fully confirm our main results, Eqs. and (|3J). 

We next turn our attention to the quantum-classical 
correspondence in phase space. We compare in 
Fig. |2] the Liouville evolution of a classical distribu- 
tion with that of the Wigner function VFi(p, q;i) — 
{2ttH)-'^ J dxexp[ipx]pi(q + x/2,q - x/2;t). The lat- 
ter is quantum-mechanically evolved from a localized 
wavepacket with the same initial location and extension 
as the classical distribution. Three quantum phase-space 
plots are shown: (i) (top right) for a free system, e = 0; 
(ii) and (iii) (bottom left and right) for a coupled system 
e = 4, in the regime 7^(i) ~ exp[— Aii]. The bottom left 
panel has a system size Ni = N2 = 512 while the bot- 
tom right panel has Ni = N2 = 2048. All plots show 
phase-space distributions after 5 kicks, a duration com- 
parable to te. Two things are clear from these figures. 
First, a coupling is necessary and sufficient to achieve 
phase-space quantum-classical correspondence. Second, 
the correspondence becomes better as we move deeper 
in the semiclassical regime, even though the interaction 
Hamiltonian vanishes in that limit ! 

One key issue is whether the observed classical entan- 
glement rate translates into a Lyapunov decoherence rate 
for systems coupled to a true environment. The latter 
differs from a coupling to a single particle in that it has 
much shorter time scales, it has a much bigger Hilbert 
space, and it cannot be initially prepared in a pure Gaus- 
sian wavepacket. We can take these conditions into ac- 
count in our semiclassical approach by considering (i) 
A2 ^ Ai, (ii) N2 Qo and (iii) taking an initial mixed 
environment density matrix ponv = X^a l</'a|^l^)(^l7 with 
(x|a) being nonoverlapping Gaussian wavepackets. The 
result is that Eq. is replaced by 

V{t) ~ aie(t > ri)exp[-Aii] -Hexp[-2rt] 

+Q{t>4^)N^\ (7) 

The dynamical Lyapunov decay of the purity seems to 
survive in the case of a particle coupled to an envi- 
ronment. We have obtained numerical confirmation of 
Eq. ITJ which we do not present here. 

We stress in conclusion that one advantage of our ap- 
proach is that V^t) is directly calculated, without the 
step of numerically integrating a differential equation for 
Pi(t). Future works should focus on multipartite entan- 
glement and decoherence by an environment consisting 
of many coupled dynamical systems. 

This work has been supported by the Swiss National 
Science Foundation. 
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FIG. 2; (Color online) Phase-space plots for a classical distribution (top left), uncoupled (top right) and coupled (bottom 
left and right, e = 4) quantum Wigner distributions, after five iterations of the kicked rotator map of Eqs. ||5J and Q. In 
all cases, the system has K\ = 3.09, and the initial distributions are Gaussian centered in the chaotic sea at {x,p) = (1,2). 
Bottom panels: Wigner functions for the quantum system coupled to a second kicked rotator with K2 = 100. One has 
2V = 13.6 > A2 ^ Ai, so that the purity behaves as V{t) ~ exp[— Ait]. The left panel has Ni = N2 = 512 and the ri ght panel 
has Ni — N2 = 2048. The presence of ghost images in the Wigner function is an artifact of the boundary conditions [lq|. 
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